Nonlinear Coherent Destruction of Tunneling 
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We study theoretically two coupled periodically-curved optical waveguides with Kerr nonlinearity. 
We find that the tunneling between the waveguides can be suppressed in a wide range of parameters 
due to nonlinearity. Such suppression of tunneling is different from the coherent destruction of tun- 
neling in a linear medium, which occurs only at the isolated degeneracy point of the quasienergies. 
We call this novel suppression nonlinear coherent destruction of tunneling. This nonlinear phe- 
nomenon can be observed readily with current experimental capability; it may also be observable 
in a different physical system, Bose-Einstein condensate. 
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Periodic driving force is an important and effective tool 
for coherently controlling quantum tunneling. This has 
been well demonstrated with a paradigmatic model, a 
free particle in a double-well potential and driven by a 
periodic external field[l[. With appropriately tuned pa- 
rameters, the periodic driving force is able not only to 
enhance tunneling [2|] -|4| but also to completely suppress 
itd-i. The latter is rather surprising and was discov- 
ered first by Grossmann et a/[5(. It is now known as 
coherent destruction of tunneling (CDT) [5| . When it oc- 
curs, a localized wave packet prepared in one well re- 
mains in the same well and does not tunnel to the other 
well. In a periodically driven system, there are Floquet 
states and associated quasienergies 0. The CDT is found 
to occur only at the isolated degeneracy point of the 
quasienergies 0, @ • 

Recently, this quantum phenomenon of CDT was ob- 
served experimen tally with two coupled periodically- 
curved waveguides [lOj (see Fig Q]). In this classical optical 
system, the Maxwellian wave mimics the quantum wave 
while the periodic driving force is achieved by bending 
the waveguides periodically. Such a waveguide system 
is an ideal laboratory system for demonstrating the co- 
herent control of quantum tunneling by periodic driving 
force. For example, tunneling enhancement has recently 
also been reported with two optical waveguides [11 1. 

In this Letter we consider a similar coupled waveguide 
system but with Kerr nonlinearity. With a well-known 
two-mode approximation, the system can be described 
by a two-mode nonlinear model with an external peri- 
odic driving force. This driving is characterized by two 
parameters, its frequency w (the inverse of the period of 
the curved waveguide) and its strength S (the curving 
magnitude of the waveguides) of the driving force. By 
numerically solving this two-mode nonlinear model, we 
find that the suppression of tunneling between the two 
coupled waveguides happens for a wide range of ratio 
S/w. This is in stark contrast to the CDT in curved lin- 
ear waveguides that occurs at an isolated point of S/w, 
where the quasienergies of the system are degenerate. 
This extension of tunneling suppression region is caused 
by nonlinearity. Therefore, we call it nonlinear coherent 



destruction of tunneling (NCDT). We find that the range 
of ratio S/w for NCDT increases steeply with nonlinear 
strength. The Floquet states and the quasienergies of 
this nonlinear model are also studied. We discover that 
there can be more than two Floquet states and quasiener- 
gies in a certain range of ratio S/w. These additional 
Floquet states form a triangle in the quasienergy levels. 
Our study reveals that these additional Floquet states 
are closely related to the NCDT. 

The current experimental capability with nonlinear 
waveguides is examined. We find that the observation 
of NCDT is well within the current experimental abil- 
ity. Note that the nonlinear two-mode model that we de- 
rived for the waveguides can also be used to describe the 
dynamics of a Bose-Einstein condensate in a double-well 
potential under a periodic modulation[l2| . This indicates 
that NCDT may also be observable with Bose-Einstein 
condensates. 




FIG. 1: Schematic drawing (not to scale) of two periodically 
curved optical waveguides placed parallel to each other. 

In a weakly guiding dielectric structure, the effective 
two-dimensional wave equation for lig ht p ropagation in 
nonlinear directional waveguides reads[13j| 
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where A is the free space wavelength of the light, xq(z) = 
Acos(2nz/A), and V(x) = [n 2 s - n 2 {x)}/(2n s ) ~ n s - 
n(x), where n(x) and n s are, respectively, the effective 
refractive index profile of the waveguides and the sub- 
strate refractive index. For the coupled waveguides as in 
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FiglU n(x) thus V(x) have a double-well structure. The 
scalar electric field is related to ip through E(x, z, t) = 
(l/2)(\n 2 \n s eoCo/2)~ 1 / 2 [tp(x, z) exp(— iuit + ikn s z) + c.c], 
where n 2 is the nonlinear refractive index of the medium, 
k = 27r/ A, lu = fcco, and cq and eo are the speed of light 
and the dielectric constant in vacuum, respectively. The 
field normalization is taken such that |'0| 2 /|ri2| gives the 
light intensity I (in W/m 2 ). By means of a Kramers- 
Henneberger transformation[3| x' = x — xq(z),z' = 
z, and 4>{x' , z') ~ ip(x' , z') exp[—i(2n s ir/ \)±o(z')x' — 

i(n s Tr/\) J* rf^ig(^)] (the dot indicates the derivative 
with respect to z'), Eq.(p} is then transformed to 
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where F(z') = n s x (z') = (A-k 2 An s / A 2 ) cos(2nz' /A) is 
the force induced by waveguide bending. It is clear that 
if we view z (or z') as time t, the above equations can be 
regarded as describing the system of a nonlinear quan- 
tum wave in a double- well potential and under a periodic 
modulation. 

We assume that the light in each waveguide of the cou- 
pler is single moded and neglect excitation of radiation 
modes. With a standard two-mode approximation [TEl 
Til . 17 1, we write 



<p(x',z') = e^ B » z ' [ci(z>i(x') + c 2 (z')u 2 (x')\ , (3) 

where u\ and u 2 are localized waves in two waveg- 
uides while the two coefficients are normalized to one, 
| C! | 2 + |c 2 | 2 = 1. E is defined as E Q — J u* 12 H ui t2 dx' . 
It is reasonable to assume that the localized wave is a 
Gaussian, u lj2 (x') = V~D exp[— (x' ± a/2) 2 /26 2 ], where 
a is the distance between the two waveguides, b is the 
half-width of each waveguide, and D is related to the 
input power of the system -P(O) as D — n 2 P(0)/(^irb). 
P(0) has the unit of W/m. The two-mode approximation 
eventually simplifies Eq.© to 



idi = -c 2 - — cos(wz')ci - x|ci| 2 ci, 
v S 



(4) 



ic 2 = -ci + — cos(wz')c 2 - x|c 2 | c 2 , (5) 

where we have set S = 8TT 3 aAn s / A 2 X, v = 
47r(J u\Hqu 2 <1x) I X, the modulation frequency w — 
27r/A, and \ = v / 27rn 2 P(0)/(A&) is an effective nonlinear 
coefficient. When 5 = 0, Eqs. id), © will be reduced 
to the well-known Jensen equation [171]. Note that -P(O) 
has the unit of W/m is because the waveguide is two 
dimensional in our theoretical model. In experiments, 
P(0) has the unit of W and the waveguides are three 
dimensional. As a result, to relate our nonlinear pa- 
rameter to the real experimental parameters, we choose 




FIG. 2: The intensity of light in the initially populated waveg- 
uide for the case of x = (dashed lines) and \/v = 0.4(solid 
lines) with (&)S/w = 1.8, (b)S/w = 2.2, (c)S/w = 2.4. Dis- 
tance z' is in units of 1/v. w/v — 10. 



X = 27m 2 P(0)/(A<7 o ff), where a c s is the effective cross- 
section of the waveguide, according to Ref. [1 81] - 

To investigate tunneling effect, we solve the two non- 
linear equations ((4]) and ([5]) numerically with the light 
initially localized in one of the two waveguides. With the 
numerical solution, we compute the intensity of the light 
staying in the initial well with P'(z') = |c*(0)ci(z') + 
c 2 (0)c 2 (z')| 2 . Three sets of our results are shown in 
Figl21[a,b,c). In the first set for S/w = 1.8, we see that 
P'(z') oscillates between zero and one for both linear case 
X = and nonlinear case x/ v = 0-4, demonstrating no 
suppression of tunneling. In the second set for S/w = 2.2, 
we see a different scenario, the oscillation of P'(z') is lim- 
ited between ~0.8 and one for the nonlinear case, showing 
suppression of tunneling, while there is no suppression 
for the linear case. In the third set for S/w — 2.4, sup- 
pression of tunneling is seen for both linear and nonlinear 
cases. Such suppression of tunneling for the linear case is 
known as coherent destruction of tunneling [5(. These nu- 
merical results demonstrate that nonlinearity can extend 
the parameter range of the suppression of tunneling. We 
call this new phenomenon nonlinear coherent destruction 
of tunneling (NCDT). 

The extension of tunneling suppression regime of ra- 
tio S/w by nonlinearity is more clearly demonstrated in 
FiglUa). In this figure, we have used localization, which 
is defined as the minimum value of P'(z'), to measure 
the suppress of tunneling. When there is large suppres- 
sion of tunneling, localization is close to one; when there 
is no suppression, localization is zero. As clearly seen in 
FigEIa), the peak of localization (solid line) for x/v = 0.4 
is much wider than the peak for xl v — 0-0 (dashed line). 
In Fig(3](b), we see the width of localization Ar increases 
almost linearly with nonlinearity x (solid line). Note that, 
analytically, CDT occurs only at isolated points. That it 
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FIG. 3: (a)Localization as a function of S/w. The solid line is 
for the nonlinear case x/ v — 0.4 and the dashed line is for the 
linear case x — 0- w/v = 10. (b) The width AF of the peak 
in (a) as a function of nonlinearity strength \/ v (solid line). 
The dashed line is for the width of the quasienergy triangle 
in FigH 



has a narrow range in FigJ3Ja) is because the evolution 
time is finite in numerical simulation. 

As is well known, the CDT is connected to the degen- 
eracy point of quasienergies in the system Q. Although 
our system is nonlinear, one can similarly define its Flo- 
quet state and quasienergy. That is, Eqs. (|4l5[) have so- 
lutions in the form of {01,02} = e~ l£Z {ci(z'), £2(2')}, 
where both c\ and £2 are periodic with period of A. These 
Floquet states and corresponding quasienergies e can be 
found numerically. We first expand the periodic functions 
61,2 in terms of Fourier series with a cutoff. After plug- 
ging them into Eqs. (|4|5| ). we obtain a set of nonlinear 
equations for the Fourier coefficients. By solving these 
equations numerically, we obtain the Floquet states and 
corresponding quasienergies e. The results are plotted in 
FigJJJ where we witness a striking difference between the 
linear and nonlinear cases. As seen in Figdta), for the 
linear case, there are two Floquet states for a given value 
of S/w and there is only one isolated degeneracy point. 
For the nonlinear case, we notice that there are four Flo- 
quet states and three quasienergies in a certain range of 
S/w with two of the Floquet states degenerate. The three 
quasienergies form a triangle in the quasienergy levels as 
seen in Figrj|b,c). Our numerical computation shows 
that the width of the quasienergy triangle increases with 
nonlinearity \ as shown in Fig[3] (dashed line). As this 
increasing trend is similar to the localization width Ar, 
this offers us the first glimpse of link between NCDT and 
the quasienergies. Since the right corner of the triangle 
can be open, we define the width of the quasienergy tri- 
angle as the horizontal distance between the left corner 
and the upper corner. 

A firm link between the NCDT and the triangle struc- 
ture in the quasienergies can be established by looking 
into the Floquet states. We focus on the Floquet states 




FIG. 4: Quasienergies at (a) \ = 0;(b) x/ v = 0.4; (c) 
X/v = 0.8. Solid lines are for numerical results obtained with 
Eqs. (|4l5p and circles for the approximation results for high 
frequencies with Eqs, (|7l8[ ). w/v = 3. 



that correspond to the lowest quasienergies in Figj4] To 
measure how the Floquet state is localized in one of the 
two waveguides, we define (|ci| 2 ) = (J Q A dz'\c\\ 2 )/k for 
a given Floquet state {01,02}- We have plotted this 
value for the lowest Floquet states in FigO In this 
figure, we see clearly that only the Floquet states on 
the quasienergy triangle are localized. This thus demon- 
strates a clear link between the quasi-energy triangle and 
the NCDT. That there are two lines in Fig[5] reflects the 
fact that there is a two-fold degeneracy for the lowest 
quasienergies on the triangle. 




FIG. 5: Intensity in the first well for every Floquet state in 
the lowest quasienergy level at x/ v — 0.4, w/v = 3. 

The triangular structure in the quasienergy is very 
similar to the energy loop discovered within the con- 
text of nonlinear Landau-Zcncr tunneling 21]. In fact, 
they are mathematically related. For high frequencies, 
w ^> max{v, x}t which is usually the case for current ex- 
periments with optical waveguides, we take advantage of 
the transformation 



Ci = c[ exp[iS , sin(wz')/2u>], 
C2 = c 2 exp[— iS sin(wz')/2w}. 



(6) 
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After averaging out the high frequency terms [Fx 
a non-driving nonlinear model, 

ic 'i = ^MS/w)c' 2 - x\c'i\ 2 c[, 



IC.-) 



-MSfw)c' x -xK\ 2 ^ 



we find 

(7) 
(8) 



where Jo is the zeroth-order Bessel function. It is clear 
from the transformation in Eq.© that the eigenstates 
of the above time-independent nonlinear equations cor- 
respond to the Floquet states of Eqs. (|4l5[ ). We have com- 
puted the eigenstates of Eqs. (|7l8[ ) and the corresponding 
eigenenergies, which are plotted as circles in FigdJ The 
consistency with the previous results is obvious. As is 
known in Ref. [2l| , the above nonlinear model admits ad- 
ditional eigenstates when \ > Jq{S/w)v. Therefore, this 
can be regarded as the condition for the extra Floquet 
states to appear for the driving nonlinear model Eqs. (|4|5[ ) 
at high frequencies. 

So far, we have focused on self-focusing materials. Our 
approach and results will be very similar if one consid- 
ers instead self-defocusing materials, for which the sign 
before the nonlinear term in Eq.JT]) should be plus. Non- 
linear coherent destruction of tunneling still occurs and 
the triangular structure also appears in the quasienergy 
levels but its direction is reversed as compared to the 
self-focusing case. 

At present the nonlinear waveguides are readily avail- 
able in labs 18j, Il9j. [20l| . We take the experimental 
parameters in Ref.[19j to estimate our theoretical val- 
ues in Eqs. (|4l5p . The wavelength of the laser light is 
A = 1.55/im, the effective cross-sectional area of the 
waveguide is a c s = 12/im 2 , the nonlinear index n,2 — 
1.2 x 10^ 13 cm 2 /W , and the shortest length for the light 
transfer from one waveguide to the other waveguide in 
the weak nonlinearity limit is L c « 2cm. With the power 
input in the waveguides P(0) ~ 100T4 7 , we have 



X = 2Trn 2 P(0)L c 

V nXcTeff 



2. 



(9) 



This shows that strong nonlinear waveguides are avail- 
able at optical labs and nonlinear coherent destruction 
of tunneling can be visualized in an optical experiment 
similar to the one in Ref.[lj|. We also want to mention 
briefly that NCDT may be applied to improve optical 
switching devices [3, Eoj . The details will be discussed 
elsewhere. 

In conclusion, we have studied the light propagation 
in a nonlinear periodically-curved waveguide directional 
coupler. We have found a new type of suppression of 
tunneling in this system, which is induced by nonlinear- 
ity and has no linear counterpart. We call it nonlinear 
coherent destruction of tunneling (NCDT) in analogy to 



a similar but different phenomenon in linear driving sys- 
tems, coherent destruction of tunneling. The NCDT oc- 
curs for an extended range of ratio S/w, where S is the 
strength of the driving and w is its frequency. We have 
also found that the NCDT is closely related to a trian- 
gular structure appeared in the quasienergy levels of the 
nonlinear system. We have also pointed out that obser- 
vation of the novel nonlinear phenomenon is well within 
the capacity of current experiments. 
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